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[ Analytical and numerical procedures are applied to show that both third and second 
order dispersion maps can be explicitly constructed and their mutual effects on the 
optical pulse propagation are analysed. In these connection it is also shown how the 
other important features such as amplification, intra-channel Raman Scattering(IRS), 
fibre loss, centre frequency of the pulse spectrum effect the propagation of pulse. Due 
to the presence of IRS, moment method is adopted which is easily reduced to usual 
Hamiltonian technique in absence of IRS. It is observed that the long time periodic 
behaviour is stabilized by the third order dispersion map. Though the introduction 
of IRS and other effects perturb the situation, one can observe that except for a very 
special case it remains stable.] 
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1 Introduction 



The study of solitary pulse propagation in optical fibre is nowadays a subject of in- 
terest not only to tlieoritical[l] and experimental physicists [2] but also to technical 
developers. Due to it's immense posibility for the world of communication, light by 
light controlled devices and ultra fast spectroscopy, various innovative ideas have been 
developed and tested. It is now a well known fact that dispersion management [3] plays 
a central role in the fabrication of long distance fibres for the optimal transmission 
of the signal. Uptill now the main emphasis is on the second order dispersion effect. 
A fundamental motivation for the dispersion mapping technique arises from the need 
to significantly suppress the phase-matched four wave mixing process, which leads 
to substantial performance penalties in wavelength division multiplexing(WDM) sys- 
tems[4]. Though dispersion management imposes a strong perturbation to the soliton 
propagation, original computational evidence of Smith et.al.[5], Gabitov et.al.[6-8] and 
Golovchenko etal.[9] suggest that solitary wave are indeed robust and can actually 
cover the long distance. It should be noted that during a given dispersion map period, 
the pulse can undergo changes in amplitude, width, and chirp. But these changes 
sampled at appropiate locations of the map can be periodic and thus, a dispersion 
managed breather is supported. 

Here we have extended the standard analysis of second order dispersion map 
to include both second and third order dispersion maps and to monitor the effect 
of amplitude, IRS, fibre loss and cetre frequency on the characteristics of the pulse 
propagation. Our method is partly analytic and partly numerical. It is observed 
from the phase plane dynamic plot that both the second and third order effect can 
mutually and effectively exist. On the other hand, we have observed that the required 
long time periodic behaviour of the pulse is altered to some extend when IRS, fibre 
loss, center frequency are taken care of. A very different phenomenon is observed 
when the second order dispersion is in the anomalous region and a dispersion map is 
considered for the third order case, which shows a drastic effect due to IRS and other 
parameters. 

Lastly it may be mentioned that since we want to include IRS and other effects 
which may not be amenable to a variational approach, we have formulated the whole 
computation by the moment method. 



2 Formulation: 

We consider the basic equation describing the pulse propagation to be the generalized 
Nonlinear Schrodinger Equation written as 

OA a iP2{z)d^A f33{z)d^A 

(1) 

where A(z,t) is the slowly varing amplitude of the pulse envelop, a accounts for 
fibre loss, is the group velocity dispersion (GVD), f]^ is the third order dispersion 
function(TOD), 7 is the nonlinear parameter responsibe for self modulation , wq is 
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the center frequency of the pulse spectrum and Tr represents the Raman parameter 
responsible for the intrapulse raman scattering(IRS). The basic idea is to treat the 
optical pulse as a particle with energy E, position T and frequency shift fl, which are 
defined via, 
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The chirp of the pulse is given as, 
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We consider the following trial wave function for the solution of Eq.(l) 
A{z, t) = \ — sech[r){t - T)] exp \ i(5{t - Tf + in{t - T) + i 
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where E(z), r]{z), P{z), il{z), T(z) and (f){z) are the soliton's energy, inverse width, 
chirp, frequency, centre and phase respectively. So, from Eqs.(2)-(6) and using Eq.(l) 
we get 
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Using the explicit form of A(z,t) given in Eq.(7) and the relations, — and 



C — Eqs. (8)-(12) can be simphfied to the following form; 
-aE 

/?3 



dE 
dz 



dT 

dz 

dn 

dz 
drj 
dz 
d£ 

dz 



+ (1 + r 



T)^' 3 



15 6ujQ 



-^Erf 



(13) 

(14) 
(15) 
(16) 
(17) 



It is now essential to specify the form of the dispersion management functions I32{z) 
and /?3(-2), which will take care of the normal and anomalous dispersion segments. 
The dispersion management function is written as; 



/?3(^)=/?2(^)=<^a + -A(C) 



(18) 



with dispersion map, 
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One can mention that it is possible to choose different form for P2 and also. The 
set of equations(13)-(17) can not be treated numerically and so we take recourse to 
Runge-Kutta-Felberge method for their solutions. 



3 Numerical Simulation: 

Before going for the complete solution we consider the case of an ideal fibre where 
Tr — 0, a — 0.0 and ^ = 0.0. In both the functions P2 and given in equation 
(18), we assume 5a — 0.0021 and Za = 1.0 where 5a is the path average dispersion and 
Za is amplifier spacing. The form of A{() is given in Fig.-l(a). After integration we 
plot the phase plane diagram in (jj, j3) plane and depict the alternative flow through 
the anomalous and normal regions by different type of lines in Figs.-l(b)-l(e) taking 
the standard form /33 = 0,/32 = I32{z) given in Eq.(18); I32{z) — I3z{z) — dispersion 
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function; (52 — —1.0, (5z{z) ^dispersion function and (52 — 1.0, (5z{z) — dispersion 
function respectively. In last two cases the shape of the phase plane are totally 
changed. In each of the figures a portion of the figures shows the reversal of the 
dynamics due to the switching of the sign of the dispersion (in fig- 1(b) and 1(c)). 
Though it is only geometric description, it shows the complete perodic character in 
absence of 1/ojq, a and Tr. Lastly we continue with the case, (52 — 0.0, (5^ =dipersion 
function, depicted in Fig-2. In Fig 2(a) we have the situation with all external effects 
put to zero and those in fig-2(b) show the closed path in presence of T^^ujq and a. 
which clearly indicates to a stronger effect of TOD (third order dispersion). 

We next analyse the explicit numerical results for the intensity(77) and chirp(/3). 
In Fig.-3 we exibit the variations of r] and (5 for the case, (52 ^dispersion function and 
(5^ = 0.0 taking different values of T/j, a and a;o(values are given in figure caption). 
The first column of Fig.-3 shows the behaviours of t] and /9 upto 10 periods of disper- 
sion map and second one exhibits for the single period after traversing nine periods of 
dispersion map. The long time periodic natures of rj and (5 are quite evident for the 
starting situation {I/coq — a — Tji — 0.0) but change slightly due to the introduction 
of these effects. We now take up the case when (52 and both are dispersion maps. 
For such a case, the variations of rj and (5 are shown in Fig.-4 and the same conclusion 
is still valid. One may note that introduction of o^o, ct, Tr have more influence than in 
the previous case. Next we consider the case (52 ~ —1 and /^a = dispersion function. 
The corresponding situations are shown in Fig.-5. Here also the starting situation 
leads to pure long time periodic nature, but remarkable changes are observed with 
the introduction of uq^u^Tr and due to this fact we have not shown the behaviour 
after nine map periods. 

At the end we depict the variations of (5 and 77 in the situation, = but (5^ — 
dispersion map. These are exibited in figures 6(a)-6(h). These figures show that even 
after a longer travelling distance (here we have considered 30 map periods), the pulse 
shows remarkable stability and the long time periodic behaviour is very much evident, 
though there may be some deviations from the required features when uiQ^a^T^ are 
introduced. This observation gives a further justification for the inclusion of third 
ordr dispersion map in the construction of fibre. 

4 Conclusion: 

Our above analysis suggests that though the basic theory of wave propagation in 
nonlinear fibre suggests the presence of higher order dispersion and nonlinearity ,for 
practical situation the effect of such term should be correctly monitored before put to 
use. In absence of other extraneous effects the closed trajectories can be observed in 
phase plane and a map constructed for both second and third order dispersion. The 
corresponding scenario of chirp and width shows that the system is highly stabilized, 
yet when IRS and other effects are taken into account and the behaviours of chirp 
and width show slight fluctuation. 
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7 Figure Caption: 



Figure-l.(a)Dispcrsion map. 

Figure-1. (77, /3) phase plane diagram with — a — Tr — 0.0 and 

b) P2— dispersion function, P3—O.O 

c) P2iz) — /33{z) = dispersion function. 

d) P2 = — 1-0,/33 = dispersion function. 

e) (32 = 1-0, /^s = dispersion function. 

Figure-2. {r),(3) phase plane diagram for the case, /?2 = 0.0, /^a = dispersion function 
with 

a) ^ = a = Tr = 0. 

h) ^ = 0.33, a = 0.0776, Tr = 1.9. 
Figure-3. The variations of 1] and /? for /?2 — dispersion function, /?3 = 0.0 with 

(a,e)7=1.0,^ = a = T^ = 0.0 

(b,f)7 = 1.0,^ = 0.2,a = T« = 0.0 

(c,g) 7 = 1.0, = 0.2, a = 0.02, Tr = 0.0 

(d,h) 7 = 1.0, = 0.2, a = 0.02, = 0.5 
Figure-4. The variations of 7] and j3 for P2 — — dispersion function with 

(a,e)7=1.0,;^ = a = T« = 0.0 

(b,f)7 = 1.0,^ = 0.2,a = T^ = 0.0 

(c,g) 7 = 1.0, = 0.2, a = 0.02, Tr = 0.0 

(d,h) 7 = 1.0, ^ = 0.2, a = 0.02, T,^ = 0.5 
Figure-5. The variations of rj and /3 for /32 = —1.0,Ps — dispersion function, with 

a) 7 = 1.0,^ = a = ri^ = 0.0 

b) 7 = 1.0,^ = 0.2,Q; = T,j = 0.0 

c) 7 = 1.0, = 0.2, a = 0.02, Tr = 0.0 

d) 7 = 1.0,J = 0.2,a = 0.02,TK = 0.5 

Figure-6. The variations of rj and /3 for (32 = 0.0 , (3^ =dispersion function with 

(a,e)7=1.0,^ = o; = r,j = 0.0 

(b,f)7 = 1.0,;j^ = 0.33,a = T^ = 0.0 

(c,g) 7 = 1.0, ^ = 0.33, a = 0.0776, Tr = 0.0 

(d,h) 7 = 1.0, ^ = 0.33, a = 0.0776, = 1.9 
(Note: The fist and second columns of Figs-3,4 correspond to the behaviours of 1] and 
(3 upto ten periods of disperison map and single period after travelling nine periods 
respectively. In first column of Fig.-6 we consider 30 map periods and second column 
shows 30th period of the map. ) 
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FIGURE - 1 




FIGURE -3 



after 9 map periods 
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FIGURE- 5 




ca 



-20 



(f) 




' ' ' ' 1 





CO. 

-1 

2 

en. 
-2 



(g 







(h) 



50 



50 



100 



100 







50 



100 



20 
10 


I 

10 

5 


I 

4 
2 


I 

2 
1 


I 

50 

■50 

I 

10 

■10 

I 

2 

-2 

I 

1 

-1 



FIGURE -6 



after 29 map periods 



JwwwwWwWWwU 



JwUwwU 



100 200 300 




100 200 300 



^VY\A/WV\AA/VWVV\AA/VW\AA/VW 



100 200 300 





100 200 300 



100 200 300 



\r\-^ 



100 200 300 



ATVAAAAAAAAAAAAAAAAAAAAAAAAAT^ 



100 200 300 




(a) 



(b) 



(c) 



(d) 



(e) 



(g) 




290 292 294 296 298 

10 I ■ ■ ■ ' 







290 292 294 296 
1.5 




290 292 294 296 298 
2 




290 292 294 296 298 
50 




290 292 294 296 298 

10 I ■ ■ ■ ' — 



(f) 



-10 



290 292 294 296 298 
0.5 




-0.5 

290 292 294 296 298 
0.5 



(h) 




-0.5 

100 200 300 290 292 294 296 



FIGURE -2 





z 



z 



